Properties of the partonic phase at RHIC 
within PHSD 



E. L. Bratkovskaya 

Institute for Theoretical Physics, University of Frankfurt, Frankfurt, Germany 
Frankfurt Institute for Advanced Study, Frankfurt, Germany 

E-mail: Elena . BratkovskayaOth . physik . uni-f rankf urt . de 

W. Cassing 

Institute for Theoretical Physics, University of Giessen, Giessen, Germany 
V. P. Konchakovski 

Institute for Theoretical Physics, University of Giessen, Giessen, Germany 
Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine 

O. Linnyk 

Institute for Theoretical Physics, University of Giessen, Giessen, Germany 
V. Ozvenchuk 

Frankfurt Institute for Advanced Study, Frankfurt, Germany 
V. Voronyuk 

Joint Institute for Nuclear Research, Dubna, Russia 
Frankfurt Institute for Advanced Study, Frankfurt, Germany 
Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine 

Abstract. The dynamics of partons, hadrons and strings in relativistic nucleus-nucleus 
collisions is analyzed within the novel Parton-Hadron-String Dynamics (PHSD) transport 
approach, which is based on a dynamical quasiparticle model for partons (DQPM) matched 
to reproduce recent lattice-QCD results - including the partonic equation of state - in 
thermodynamic equilibrium. The transition from partonic to hadronic degrees of freedom is 
described by covariant transition rates for the fusion of quark-antiquark pairs or three quarks 
(antiquarks), respectively, obeying flavor current-conservation, color neutrality as well as energy- 
momentum conservation. In order to explore the space-time regions of 'partonic matter' the 
PHSD approach is applied to nucleus-nucleus collisions from SPS to RHIC energies. Detailed 
comparisons are presented for hadronic rapidity spectra and transverse mass distributions. The 
traces of partonic interactions are found in particular in the elliptic flow of hadrons as well as 
in an approximate quark-number scaling at the top RHIC energy. 



1. Introduction 

The 'Big Bang' scenario implies that in the first micro-seconds of the universe the entire state has 
emerged from a partonic system of quarks, antiquarks and gluons - a quark-gluon plasma (QGP) 
- to color neutral hadronic matter consisting of interacting hadronic states (and resonances) 
in which the partonic degrees of freedom are confined. The nature of confinement and the 
dynamics of this phase transition has motivated a large community for several decades and is 
still an outstanding question of todays physics. Early concepts of the QGP were guided by the 
idea of a weakly interacting system of partons which might be described by perturbative QCD 
(pQCD). However, experimental observations at the Relativistic Heavy Ion Collider (RHIC) 
indicated that the new medium created in ultrarelativistic Au+Au collisions is interacting more 
strongly than hadronic matter and consequently this concept had to be severely questioned. 
Moreover, in line with theoretical studies in Refs. [H El E] the medium showed phenomena of 
an almost perfect liquid of partons [H [5] as extracted from the strong radial expansion and the 
scaling of elliptic flow V2(pt) of mesons and baryons with the number of constituent quarks and 
antiquarks [I]. 

The question about the properties of this (nonperturbative) QGP liquid is discussed 
controversially in the literature and dynamical concepts describing the formation of color neutral 
hadrons from colored partons are scarce. A fundamental issue for hadronization models is the 
conservation of 4- momentum as well as the entropy problem because by fusion/coalescence of 
massless (or low constituent mass) partons to color neutral bound states of low invariant mass 
(e.g. pions) the number of degrees of freedom and thus the total entropy is reduced in the 
hadronization process [3 [8] . This problem - a violation of the second law of thermodynamics 
as well as the conservation of four-momentum and flavor currents - has been addressed in Ref. [9 J 
on the basis of the DQPM employing covariant transition rates for the fusion of 'massive' quarks 
and antiquarks to color neutral hadronic resonances or strings. In fact, the dynamical studies 
for an expanding partonic fireball in Ref. [9j suggest that the latter problems have come to a 
practical solution. 

A consistent dynamical approach - valid also for strongly interacting systems - can be 
formulated on the basis of Kadanoff-Baym (KB) equations [10] or off-shell transport equations 
in phase-space representation, respectively [10J [TTJ [12]. In the KB theory the field quanta are 
described in terms of dressed propagators with complex selfenergies. Whereas the real part of 
the selfenergies can be related to mean-field potentials (of Lorentz scalar, vector or tensor type), 
the imaginary parts provide information about the lifetime and/or reaction rates of time- like 
'particles' [13j . Once the proper (complex) selfenergies of the degrees of freedom are known the 
time evolution of the system is fully governed by off-shell transport equations (as described in 
Refs. [10\ I13|). The determination/extraction of complex selfenergies for the partonic degrees 
of freedom has been performed before in Refs. [HI Q2] by fitting lattice QCD (1QCD) 'data' 
within the Dynamical QuasiParticle Model (DQPM). In fact, the DQPM allows for a simple 
and transparent interpretation of lattice QCD results for thermodynamic quantities as well as 
correlators and leads to effective strongly interacting partonic quasiparticles with broad spectral 
functions. For a review on off-shell transport theory and results from the DQPM in comparison 
to 1QCD we refer the reader to Ref. |13j . 

The actual implementations in the PHSD transport approach have been presented in detail 
in Refs. |16t 117]. Here we report again on the actual description of hadronization (Section [2]). 
We also present results for rapidity distributions, transverse mass spectra and elliptic flow for 
heavy ion collisions at SPS (Section [3j) and RHIC energies (Section 2|) in comparison to data 
from the experimental collaborations. 



2. Hadronization in PHSD 

The hadronization, i.e. the transition from partonic to hadronic degrees of freedom, is described 
in PHSD by local covariant transition rates as introduced in Ref. [9] e.g. for q + q fusion to a 
meson m of four-momentum p = (w,p) at space-time point x = (t, x): 

dN m (x,p) 4 ^z4f x q + x > 



Tr q Tr q 5 4 (p - p q - p q ) <5 4 ~^ Xq - x^j 



d 4 xSp 

XUJ q Pq{Pq) ^q Pq{Pq) \Vqqf W m (x q - Xq, (p q - Pq)/2) 

xN q (xq,p q ) N q (xq,p q ) ^(flavor, color). (1) 
In Eq. dU) we have introduced the shorthand notation, 

where J2j denotes a summation over discrete quantum numbers (spin, flavor, color); Nj(x,p) is 
the phase-space density of parton j at space-time position x and four-momentum p. In Eq. ([TJ 
5 (flavor, color) stands symbolically for the conservation of flavor quantum numbers as well as 
color neutrality of the formed hadron m which can be viewed as a color-dipole or 'pre-hadron'. 
Furthermore, v q q-(p p ) is the effective quark-antiquark interaction from the DQPM (displayed in 
Fig. 10 of Ref. [15]) as a function of the local parton (q + q + g) density p p (or energy density). 
Furthermore, W m (x,p) is the dimensionless phase-space distribution of the formed 'pre-hadron', 
i.e. 

W m (Z, P!; ) = exp exp (2b 2 (pj - (M q - M q ) 2 /4)) (3) 

with £ = x\ — X2 = x q — x q and p^ = (pi — p2)/2 = (p q — Pq)/2- The width parameter b is 
fixed by \J (r 2 ) = b = 0.66 fm (in the rest frame) which corresponds to an average rms radius 
of mesons. We note that the expression ([3]) corresponds to the limit of independent harmonic 
oscillator states and that the final hadron-formation rates are approximately independent of the 
parameter b within reasonable variations. By construction the quantity Q is Lorentz invariant; 
in the limit of instantaneous 'hadron formation', i.e. £° = 0, it provides a Gaussian dropping 
in the relative distance squared (ri — Y2) 2 ■ The four-momentum dependence reads explicitly 
(except for a factor 1/2) 

(Et - E 2 f - ( Pl - P2 ) 2 - (Mi - M 2 f < (4) 

and leads to a negative argument of the second exponential in ([3]) favoring the fusion of partons 
with low relative momenta p q — pq = p\ — p 2 - 

Related transition rates (to Eq. ([TJ) are defined for the fusion of three off-shell quarks 
(qi + Q2 + Q3 ^ B) to color neutral baryonic (B or B) resonances of finite width (or strings) 
fulfilling energy and momentum conservation as well as flavor current conservation using Jacobi 
coordinates (cf. Ref. [T6]). 

On the hadronic side the PHSD transport approach includes explicitly the baryon octet and 
decouplet, the 0~- and l~-meson nonets as well as selected higher resonances as in HSD |18j . 
Hadrons of higher masses (> 1.5 GeV in case of baryons and > 1.3 GeV in case of mesons) are 
treated as 'strings' (color-dipoles) that decay to the known (low-mass) hadrons according to the 
JETSET algorithm [19] . We discard an explicit recapitulation of the string decay and refer the 
reader to the original work [19] or Ref. |20j . 



3. Application to nucleus-nucleus collisions at SPS energies 

In this Section we employ the PHSD approach to nucleus-nucleus collisions at moderate 
relativistic energies. It is of interest, how the PHSD approach compares to the HSD [18] model 
(without explicit partonic degrees-of- freedom) as well as to experimental data. In Fig. [1] we show 
the transverse mass spectra of tt~ , K + and K~ mesons for 7% central Pb+Pb collisions at 40 
and 80 A-GeV and 5% central collisions at 158 A-GeV in comparison to the data of the NA49 
Collaboration |21j . Here the slope of the tt~ spectra is only slightly enhanced in PHSD relative 
to HSD which demonstrates that the pion transverse motion shows no sizeable sensitivity to 
the partonic phase. However, the transverse mass spectra are substantially hardened with 
respect to the HSD calculations at all bombarding energies - i.e. PHSD is more in line with 
the data - and thus suggests that partonic effects are better visible in the strangeness-degrees 
of freedom. The hardening of the kaon spectra can be traced back to parton-parton scattering 
as well as a larger collective acceleration of the partons in the transverse direction due to the 
presence of repulsive vector fields for the partons. The enhancement of the spectral slope for 
kaons and antikaons in PHSD due to collective partonic flow shows up much clearer for the 
kaons due to their significantly larger mass (relative to pions). We recall that in Refs. [22J 
the underestimation of the slope by HSD (and also UrQMD) had been suggested to be a 
signature for missing partonic degrees of freedom; the present PHSD calculations support this 
early suggestion. Moreover, the PHSD calculations for RHIC energies show a very similar trend 
- the inverse slope increases by including the partonic phase. 




Figure 1. The 7r~, K + and K~ transverse mass spectra for central Pb+Pb collisions at 40, 
80 and 158 A-GeV from PHSD (thick solid lines) in comparison to the distributions from HSD 
(thin solid lines) and the experimental data from the NA49 Collaboration |21j . 



The strange antibaryon sector is of further interest since here the HSD calculations have 
always underestimated the yield [23J. Our detailed studies in Ref. [16] show that the HSD 
and PHSD calculations both give a reasonable description of the A + S° yield of the NA49 
Collaboration [24] : both models underestimate the NA57 data |25j by about 30%. An even larger 
discrepancy in the data from the NA49 and NA57 Collaborations is seen for (A + E°) / 'N woun d) 
here the PHSD calculations give results which are in between the NA49 data and the NA57 data 
whereas HSD underestimates the (A + S°) midrapidity yield at all centralities. 

The latter result suggests that the partonic phase does not show up explicitly in an enhanced 
production of strangeness (or in particular strange mesons and baryons) but leads to a different 
redistribution of antistrange quarks between mesons and antibaryons. In fact, as demonstrated 



in Ref. |16j . we find no sizeable differences in the double strange baryons from HSD and PHSD 
- in a good agreement with the NA49 data - but observe a large enhancement in the double 
strange antibaryons for PHSD relative to HSD. 



4. Application to nucleus-nucleus collisions at RHIC energies 

In this section we continue the comparison of the PHSD transport approach to the experimental 
data from the RHIC collaborations as well as to the correspondent HSD results [T7]. We find 
the rapidity distributions of the charged mesons to be slightly narrower in PHSD than those 
from HSD and actually closer to the experimental data. Also note that there is slightly more 
production of mesons in PHSD than in HSD while the number of charged pions is slightly 
lower. The actual deviations between the PHSD and HSD spectra are not dramatic but more 
clearly visible than at SPS energies (cf. Ref. |16|). Nevertheless, it becomes clear from Fig. [2]that 
the energy transfer in the nucleus-nucleus collision from initial nucleons to produced hadrons - 
reflected dominantly in the light meson spectra - is rather accurately described by PHSD. Fig. [2] 
also demonstrates that the longitudinal motion is well understood within the PHSD approach. 
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Figure 2. The rapidity distribution of ir + (upper part, l.h.s.), K + (lower part, l.h.s.), tt~ (upper 
part, r.h.s.) and K~ (lower part, r.h.s.) for 5% central Au+Au collisions at yfs = 200 GeV from 
PHSD (solid red lines) in comparison to the distribution from HSD (dashed blue lines) and the 
experimental data from the RHIC Collaborations |26 [ 127 } [28] . 



Independent information on the active degrees of freedom is provided by transverse mass 
spectra of the hadrons especially in central collisions. The actual results for RHIC energies 
are displayed in Fig. [3] where we show the transverse mass spectra of 7r~, K + and K~ mesons 
for 5% central Au+Au collisions at y/s = 200 GeV in comparison to the data of the RHIC 
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Figure 3. The tt , K + and K transverse mass spectra for 5% central Au+Au collisions 
at y/s = 200 GeV from PHSD (solid red lines) in comparison to the distributions from 
HSD (dashed blue lines) and the experimental data from the BRAHMS, PHENIX and STAR 
Collaborations [261 EH EE] at midrapidity. 

Collaborations [261 EH EE] • Here the slope of the ir~ spectra is slightly enhanced in PHSD (solid 
red lines) relative to HSD (dashed blue lines) which demonstrates that the pion transverse mass 
spectra also show some sensitivity to the partonic phase (contrary to the SPS energy regime). 
The transverse mass spectra are substantially hardened with respect to the HSD calculations 
- i.e. PHSD is more in line with the data - and thus suggest that partonic effects are better 
visible in the strangeness degrees-of-freedom. The hardening of the kaon spectra can be traced 
back to parton-parton scattering as well as a larger collective acceleration of the partons in 
the transverse direction due to the presence of the repulsive mean-field for the partons. The 
enhancement of the spectral slopes for kaons and antikaons in PHSD due to collective partonic 
flow shows up much clearer for the kaons due to their significantly larger mass (relative to pions). 

Of additional interest are the collective properties of the strongly interacting system which 
are explored experimentally via the elliptic flow 



of hadrons as a function of centrality, rapidity y, transverse momentum or transverse kinetic 
energy per participating quarks and antiquarks. We note that the reaction plane in PHSD is 
given by the x — z plane with the z-axis in beam direction. 

We start in Fig. U] with the elliptic flow V2 (for Au+Au collisions at the top RHIC energy) 
as a function of the centrality of the reaction measured by the number of participating nucleons 
Npart- The solid (red) line stands for the results from PHSD which is compared to the data 
for charged particles from the PHOBOS Collaboration |30j. The dashed blue line refers to 
the corresponding results for V2 from HSD (taken from Ref. |29j). The momentum integrated 
results in the pseudo-rapidity window \r)\ < 1 from PHSD compare well to the data from Ref. [30] 
whereas the HSD results clearly underestimate the elliptic flow as pointed out before in Ref. [29] . 
The relative enhancement of vi in PHSD with respect to HSD can be traced back to the high 
interaction rate in the partonic phase and to the repulsive mean field for partons. We note in 
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Figure 4. The elliptic flow V2 for Au+Au collisions at the top RHIC energy y/s = 200 GeV as 
a function of the centrality measured by the number of participating nucleons N par t- The solid 
(red) line stands for the results from PHSD whereas the dashed (blue) line represents the results 
from HSD (from Ref. [29]). The data are taken from the PHOBOS Collaboration [30] and 
correspond to momentum integrated events in the pseudo-rapidity window |r/| < 1 for charged 
particles. The shaded band signals the statistical uncertainties of the PHSD calculations. 




Figure 5. Left: The hadron elliptic flow V2 for inclusive Au+Au collisions as a function of 
the transverse momentum pj- (in GeV) for different hadrons in comparison to the data from the 
STAR |31j and PHENIX Collaborations [32] within the same rapidity cuts. Right: The elliptic 
flow i>2 - scaled by the number of constituent quarks n q - versus the transverse kinetic energy 
divided by n q for different hadron species in comparison to the data from the STAR [31] and 
PHENIX [32] Collaborations. 

passing that PHSD calculations without mean fields only give a tiny enhancement for the elliptic 
flow relative to HSD. 

Fig. [5] (l.h.s.) shows the final hadron V2 versus the transverse momentum px for different 




particle species in comparison to the data from the STAR [31J and PHENIX Collaborations |32j. 
We observe a mass separation in pr as well as a separation in mesons and baryons for > 
2 GeV roughly in line with data. The elliptic flow of mesons is slightly underestimated for 
Pt > 2 GeV in PHSD which is opposite to ideal hydrodynamics which overestimates V2 at high 
transverse momenta. On the other hand, the proton (and antiproton) elliptic flow is slightly 
overestimated at low transverse momenta pt < 1-5 GeV. 

A further test of the PHSD hadronization approach is provided by the 'constituent quark 
number scaling' of the elliptic flow V2 which has been observed experimentally in central Au+Au 
collisions at RHIC [H \'6'2\ [55]. In this respect we plot V2/n q versus the transverse kinetic energy 
per constituent parton, 

KE T = , (6) 

n q 

with niT and m denoting the transverse mass and actual hadron mass, respectively. For 
mesons we have n q = 2 and for baryons/antibaryons n q = 3. The results for the scaled 
elliptic flow are shown in Fig. [5] (r.h.s.) in comparison to the data from the STAR [31] and 
PHENIX Collaborations [32J for different hadrons and suggest an approximate scaling. For 
KEt > 0.5 GeV there is a tendency to underestimate the experimental measurements for 
A, S, A, £ baryons which we attribute to an underestimation of interaction terms in PHSD for 
high momentum hadrons. In this respect we recall that the momentum independence of the 
quasiparticle width 7 and mass M is presently a rough approximation and has to be refined. 
Due to the limited statistics especially in the baryonic sector with increasing pr this issue will 
have to be re-addressed with high statistics in future. 



5. Summary 

In this contribution we have addressed relativistic collisions of Pb+Pb at SPS energies and 
Au+Au collisions at RHIC energies in the PHSD approach which includes explicit partonic 
degrees of freedom as well as dynamical local transition rates from partons to hadrons (pQ) . The 
hadronization process conserves four-momentum and all flavor currents and slightly increases the 
total entropy since the 'fusion' of rather massive partons dominantly leads to the formation of 
color neutral strings or resonances that decay microcanonically to lower mass hadrons. Since this 
dynamical hadronization process increases the total entropy the second law of thermodynamics 
is not violated (as is the case for simple coalescence models incorporating massless partons). 

The PHSD approach has been also applied to nucleus-nucleus collisions from 40 to 160 A-GeV 
as well as for RHIC energies in order to explore the space-time regions of 'partonic matter' |16j . 
We have found that even central collisions at the top SPS energy of ~158 A- GeV show a large 
fraction of non-partonic, i.e. hadronic or string-like matter, which can be viewed as a 'hadronic 
corona' [3l]- This finding implies that neither purely hadronic nor purely partonic 'models' 
can be employed to extract physical conclusions in comparing model results with data. On the 
other hand - studying in detail Pb+Pb reactions at SPS energies in comparison to the data [16] 
- it is found that the partonic phase has only a very low impact on the longitudinal rapidity 
distributions of hadrons but a sizeable influence on the transverse-mass distribution of final 
kaons due to the partonic interactions. The most pronounced effect is seen on the production of 
multi-strange antibaryons due to a slightly enhanced ss pair production in the partonic phase 
from massive time-like gluon decay and a more abundant formation of strange antibaryons in 
the hadronization process. This enhanced formation of strange antibaryons in central Pb+Pb 
collisions at SPS energies by hadronization supports the early suggestion by Braun-Munzinger 
and Stachel |35[ [36] in the statistical hadronization model - which describes well particle ratios 
from AGS to RHIC energies. 

At RHIC energies the PHSD calculations show also a good reproduction of the hadron 
transverse mass and rapidity spectra. Furthermore, the elliptic flow vi is well described for 



Au+Au reactions at y/s = 200 GeV as a function of centrality as well as of transverse momenta 
up to pt — 1.5 GeV/c. Due to the local transition rates from partons to hadrons (pQ) the PHSD 
approach also gives approximate quark number scaling of the elliptic flow as found experimentally 
by the RHIC Collaborations. 
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